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Рассмотрена статическая упругопластическая задача для двух контактирующих с натягом неоднородных по физико-механичес- 
ким характеристикам цилиндрических тел. Для вычисления контактного давления предложен алгоритм, который реализует 
идею определения положения контактной поверхности путем итерационного выравнивания одинаковых в положении равнове- 
сия радиальных напряжений контактирующих тел. Численное решение осуществляется вариационно-разностным методом с ис- 
пользованием соотношений деформационной теории пластичности и уравнений Коши. Получены параметры напряженно-де- 
формированного состояния обоих тел, а также распределение реакции ведущего устройства по его ширине для разных матери- 
алов. Результаты численного анализа приведены в виде графиков. 


1. Введение 

Решение динамической задачи определения 
напряженно-деформированного состояния (НДС) 
метаемого элемента (МЭ) в трубе метательной ус- 
тановки (МУ) можно проводить в квазистатичес- 
кой постановке, обоснование которой представле- 
но в [1-3]. При этом на него действуют давление 
газов, реакция ведущего устройства (ВУ) и появля- 
ющиеся осевые массовые силы, которые являются 
следствием использования связанной с МЭ неи- 
нерциальной системы координат г, ср, і. 

Среди перечисленных силовых нагрузок, 
действующих на МЭ при метании, реакция ВУ име- 
ет значительную величину и играет важную роль. От 
ее воздействия в корпусе МЭ могут возникать оста- 
точные пластические деформации. Поэтому знание 
величины реакции ВУ необходимо как при констру- 
ировании МЭ, так и для проведения прочностных 
расчетов корпуса МЭ и трубы МУ. Для определения 
реакции ВУ применяются теоретико-эксперимен- 
тальные методики, использующие результаты экспе- 
римента и аналитические методы, основанные на 
соотношениях сопротивления материалов и (или) 
уравнениях теории оболочек. Применяются также 
методики, которые с помощью реализующих урав- 
нения теории упругости и пластичности современ- 
ных численных методов в сочетании с результатами 
эксперимента позволяют не только определять иско- 
мую реакцию ВУ, но и проводить анализ НДС МЭ. 

2. Алгоритм определения реакции ВУ 

Представленный алгоритм определения реакции 
ВУ р изложен для совместного осесимметричного 
упругопластического деформирования системы 
"труба МУ - ВУ - корпус МЭ - снаряжение" (рис. 1 ). 

Суть реализуемой задачи взаимодействия трубы 
и корпуса МЭ формулируется следующим образом: 
найти расположение контактной поверхности АВ 
(рис. 1, б) этой системы при посадке их друг на дру- 
га с натягом, равным величине 25. В модельной осе- 
симметричной задаче рассматриваются участки тру- 
бы МУ и корпуса МЭ, прилегающие к ВУ. На рис. 1, 
й эти элементы представлены до их посадки друг на 
друга. МЭ (верхний цилиндр на рис. 1, а) состоит из 


корпуса (1), снаряжения (2) и ВУ (3). Нижний ци- 
линдр (4) представляет трубу МУ. При рассмотре- 
нии такой модели взаимодействия трубы с корпусом 
МЭ через ВУ принимаются следующие допущения: 
1) материал ВУ в процессе посадки не уносится с его 
поверхности; 2) внутренняя поверхность трубы яв- 
ляется цилиндрической; 3) не учитывается давление 
газа и действие массовых сил на трубу и корпус МЭ 
ввиду малости этих нагрузок по сравнению с реак- 
цией р на данном этапе движения МЭ. 

Оценочные расчеты показывают, что при дли- 
тельности процесса метания 0,06 с (что приемлемо 
для пороховых систем) и скорости упругого возму- 
щения в металле МЭ порядка 5000 м/с время про- 
бега волны по длине МЭ (это максимум 0,5 м) сос- 
тавит около 0,0001 с, а за время метания получится 
600 пробегов. В этих условиях нет смысла рассмат- 
ривать волновую (динамическую) постановку, и, 
учитывая все вышесказанное, задачу можно решать 
в статической постановке. 
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Рис. 1 . Расчетная схема конструкции до (а) и после (6) по- 
садки друг на друга 
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Решение задачи о посадке с натягом для кон- 
тактирующих по всей длине боковой поверхности 
цилиндров в упругой постановке получено в виде 
конечных аналитических зависимостей в курсах 
"Сопротивления материалов" и "Теории упругости" 
[2, 4] при анализе составных цилиндров. При этом 
величина натяга А (разность диаметров посадоч- 
ных поверхностей наружного и внутреннего цили- 
ндров) в цилиндрической системе координат г, і 
удовлетворяет следующему тождеству: 

А = 2 (щ-щ), (1) 

где щ, щ — радиальные перемещения посадочных 
поверхностей внутреннего и наружного цилиндров 
соответственно; щ < 0, а щ > 0. На посадочной (кон- 
тактной) поверхности АВ, определяемой переме- 
щениями щ, щ, в условиях равновесия имеет место 
равенство радиальных напряжений 

а л - а й = 0, (2) 

где а А , о а - сжимающие напряжения на посадоч- 
ных поверхностях внутреннего и наружного цили- 
ндров соответственно. 

Для решения задачи предлагается алгоритм, 
протестированный на решении упругой задачи о 
посадке двух цилиндров по всей их длине, для ко- 
торой имеется аналитическое решение, и апроби- 
рованный на решении упругой задачи о посадке с 
натягом двух цилиндров для случая контакта меж- 
ду ними по части поверхности [5, 6]. На исследуе- 
мую область наносится конечно-разностная сетка. 
Задавая в /узлах контактной поверхности ВУ ради- 
альные перемещения и и (і =1, 2, ..., /), определяем 
НДС МЭ, в том числе радиальные напряжения а ЛІ 
на внешней контактной поверхности АВ ВУ. Далее, 
используя уравнение (1), находим величины пере- 
мещений и ѵ на контактной поверхности трубы МУ 
и гі = А/2 + и и , і =1,2, ...,/, (3) 

и, определяя ее НДС, вычисляем радиальные напря- 
жения а т на внутренней контактной поверхности АВ 
трубы МУ. Затем проверяется условие (2). Если оно не 
выполняется с заранее заданной точностью, прово- 
дится корректировка заданных на контактной пове- 
рхности перемещений, допустим и и , по формуле 

иГ=и<‘>(1-сф сгЦ<тЯ-(&)/(о®+о®), і=1,2,...,І, (4) 
где к - индекс, определяющий номер итерации; 
с, - корректирующие добавки, на величину которых 
изменяются радиальные перемещения щ;, В - мате- 
матическая вязкость, вводимая для ускорения схо- 
димости. Величина с, положительна, если |ст г1 ,|>|ст й ,|, 
и отрицательна - в противоположном случае. 

Корректирующая формула (4) работает следую- 
щим образом. Если в /-ом узле контактной поверх- 
ности вычисленные радиальные напряжения удов- 
летворяют неравенству |<т^|>|ст^ ) )|, использование 
формулы (4) приводит к уменьшению абсолютной 
величины перемещения и( Ж) и, согласно выраже- 
нию (3), увеличению перемещения и$ + '\ В резуль- 
тате расчета НДС двух тел с новыми значениями 
заданных на контактной поверхности перемеще- 
ний и‘Е'\ г4* +1) абсолютная величина напряжения 


<т,[‘ +1) несколько уменьшится, напряжения о$ +1) - 
увеличится, а разница А(т(‘ +1) =|ст,!- +1) -о-,)- +1) | станет 
ближе к нулю. Аналогично происходит выравнива- 
ние напряжений в случае, когда |ст г )- ) |<|о;)‘ ) |. Итак, с 
помощью формулы (4) корректируются, по сути 
дела, все 2/ радиальные перемещения контактных 
поверхностей обоих тел. 

Описанный процесс определения НДС каждого 
из двух контактирующих тел повторяется до тех 
пор, пока условие (2) не будет выполнено. Итогом 
этого итерационного процесса является искомая 
величина реакции ВУ /»<т н ~сг й . Представленный 
алгоритм аналогичен предложенному в [7] подходу. 

Изложенная методика в части определения упру- 
гопластического осесимметричного НДС обоих ци- 
линдрических тел реализовывалась с помощью вари- 
ационно-разностного метода (ВРМ) с использовани- 
ем соотношений деформационной теории пластич- 
ности и уравнений Коши [8-1 1]. При дискретизации 
функционала энергии для производных применялись 
так называемые естественные аппроксимации [12]. 

3. Численный пример 

Представленный алгоритм был применен для 
решения модельной задачи о взаимодействии ВУ 
МЭ и трубы МУ. Параметры, описывающие геомет- 
рию трубы и части корпуса МЭ, принимались сле- 
дующими: 5 = 0,025 см; внутренний радиус трубы 
7?! =7,775 см; внешний радиус трубы К 2 =10,5 см; ра- 
диус корпуса МЭ В =7,5 см; радиус МЭ по ВУ 
Т ?4 =К { +5 =7,8 см; длина / =10,0 см, ширина ВУ 
Ц =\, 4 см; / 2 =2,0 см. Параметр / был определен пу- 
тем численного эксперимента. В качестве / приня- 
та такая его величина, увеличение которой уже не 
влияет на результаты расчетов. Физико-механичес- 
кие характеристики (ФМХ) материала трубы МУ, 
корпуса МЭ, ВУ и снаряжения приведены в табли- 
це, где Е - модуль упругости; Е л - модуль упрочне- 
ния; /и - коэффициент Пуассона; е, - деформация 
начала текучести. Материалом трубы МУ является 
сталь с довольно высоким пределом текучести. 

Таблица. Физико-механические характеристики конструкции 


ФМХ 

Труба 

МУ 

Корпус 

МЭ 

Снаря- 

жение 

ву ! 

Медь 

Пластик 

Е, МПа 

2,02х10 5 

2,02х10 5 

875 

10 3 

8,16x1 0 4 

Е,, МПа 

8,0х10 3 

8,0х10 3 

250 

8,0х10 3 

2,04х10 3 

Р 

0,333 

0,333 

0,350 

0,330 

0,230 

е % 

0,00700 

0,00173 

0,00400 

0,00173 

0,01500 


Задача решалась при следующих граничных ус- 
ловиях для обоих торцов трубы МУ и МЭ: 

1=0, Е. и(г,0)=и(г,Е)=\ѵ(г,0)=\ѵ(г,0)=0, (5) 

1=0, Е. \ѵ(г,0)=\ѵ(г,Е)=0, о г ,(г,0)=о г {г,Е)=0, (6) 

т.е. жесткой заделки и возможности перемещения их 
в радиальном направлении соответственно, ибо ре- 
альные условия на границе расчетной области таки- 
ми двумя постановками (5) и (6) берутся "в вилку". 
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На оси вращения От задаются условия симметрии: 
г = 0: и(0,і) = 0 , сг^(0,2) = 0 . 

Внешняя поверхность трубы свободна от нагрузок: 

<5 г (ЯъІ) = 0, <т,.( К ъ у) = 0, 

также, как и ее внутренняя поверхность г=К } и 
внешняя поверхность корпуса МЭ г=К ь кроме 
участков длиной Ь ь по которым осуществляется 
контакт трубы МУ и МЭ. На этих участках реализу- 
ются следующие условия: 

{1+Ц)/2>&(1-Ц)/2: и 2 (г ь і)-щ(г к ,д=8, 

<? г] (г к ,і)-о й {г к ,і)=Ъ, а Г7] (г ь і)=а пг (г к ,і)={), 
где о пі { г ъі)> & г а( г кЛ) ~ сдвиговые напряжения на 
контактной поверхности ВУ и трубы соответствен- 
но; г к - искомые радиальные координаты узлов кон- 
тактной поверхности. На контактных поверхностях 
корпуса МЭ со снаряжением и ведущим устрой- 
ством задаются условия непрерывного контакта. 

Приведенные ниже результаты численного ана- 
лиза получены для варианта (5) граничных условий 
на торцах 2=0, 1 и, ввиду симметрии НДС рассмат- 
риваемой системы относительно координаты 
1 =Ь/ 2, представлены лишь для половины длины 
конструкции. ВУ изготовлено из меди. 



Рис. 2. Распределение напряжений у внутренней поверхнос- 
ти трубы (г = 7,8 см) 

На рис. 2 изображены зависимости радиальных 
ст г , осевых < 7 ,, окружных ег напряжений и интенсив- 
ностей напряжений <т, в трубе от осевой координа- 
ты і для г = 7,8 см, т.е. у внутренней поверхности 
трубы МУ. О характере распределения и величине 
искомой реакции р можно судить по поведению 
кривой радиального напряжения о г на контактной 
поверхности для Т/2>г>(Т-Т,)/2. Величина напря- 
жения < 7 „ почти нулевая на отдалении от места кон- 
такта с ВУ, приближаясь к контактной поверхнос- 
ти, достигает наибольшего значения в точках, соот- 
ветствующих координатам <%(/.-/.,)/2, и несколько 


уменьшается к средней по ширине части %=Ь/2 ВУ. 
Такое распределение контактного давления р к ха- 
рактерно для случая посадки с натягом друг на дру- 
га двух цилиндров для случая контакта между ними 
по части поверхности [2, 4]. Осредненная по шири- 
не ВУ величина реакции рРю - 540 МПа. Распределе- 
ние напряжений в корпусе МЭ для цилиндричес- 
кой поверхности с радиусом г= 6,9 см, в том числе 
под канавкой для ВУ, представлено на рис. 3. 



Рис. 3. Распределение напряжений в корпусе МЭ для пове- 
рхности г = 6,9 см 


На рис. 4 приведены эпюры напряжений по тол- 
щине МЭ и трубы МУ в их средней по длине части 
(%=Ь/ 2). Заштрихованные зоны обозначают снаря- 
жение (4,5 см>г>0) и ВУ (7,8 см>г>7,1 см), между 
которыми находится корпус (не заштрихованная 
часть). Радиальные напряжения а г на контактной 
поверхности в МЭ и трубе достигают величины око- 
ло -540 МПа. Такую же величину имеет напряжение 
<у г у внутренней поверхности трубы для г~5 см на 
рис. 2. В снаряжении напряжения малы. При пере- 
ходе через разделяющую снаряжение и корпус МЭ 
поверхность все напряжения, кроме 8 Г , возрастают 
примерно на два порядка и больше. Так, величина 
интенсивности напряжений а; с 2 МПа в снаряже- 
нии возрастает примерно до 350 МПа в корпусе МЭ. 



Рис. 4. Распределение напряжений по толщине корпуса МЭ 
и трубы МУ для і = 1/2. Слева от оси вращения сис- 
темы на право изображены последовательно снаря- 
жение, корпус МЭ, ВУ и труба МУ. Снаряжение и ВУ 
показаны штриховкой 
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Анализ полученных напряжений для трубы поз- 
волил сделать вывод об увеличении их значений 
при приближении к центру конструкции для і=Ь / 2 
и уменьшению радиальных напряжений а; при 
движении от внутренней поверхности трубы, где 
имеет место контакт с ВУ, к свободной от нагрузок 
наружной поверхности, что согласуется с реальной 
картиной распределения напряжений. 

Радиальные перемещения внутреннего контура 
трубы МУ (кривая 1) и внешнего контура МЭ (кри- 
вая 2) изображены на рис. 5. Их значения для 
і/2>2>(Х-і!)/ 2 удовлетворяют тождеству (1). Ана- 
лиз результатов позволил выделить зоны пласти- 
ческого деформирования конструкции: всё ВУ 
(е,тах = 0,02) и часть корпуса МЭ под ВУ. 

16 

12 

8 

4 


0 

-4 

-8 

-12 


Рис. 5. Радиальные перемещения: 1) внутреннего контура 
трубы МУ; 2) внешнего контура корпуса МЭ 

На рис. 6 представлены распределения реакции р 
по ширине ВУ для разных материалов. Зависимость 
р(ф), описанная кривой 2, получена для медного ВУ. 
Кривая 1 получена также для медного ВУ, но при от- 
сутствии снаряжения. Величина реакции в этом слу- 
чае получается несколько большей, нежели при на- 
личии снаряжения. Для случая, когда ВУ выполнено 
из пластика на основе полиэтилентерефталата 
(см. таблицу), распределение реакции р описывает- 
ся кривой 3. Кривой 4 представлено распределение 
реакции для ортотропного материала ВУ. 

Осевые перемещения в узлах на контактной по- 
верхности МУ и ВУ имеют величину одного поряд- 
ка, и на один-два порядка меньше (по абсолютной 
величине), нежели радиальные перемещения. Это 
соотношение между величинами перемещений име- 
ет место для таких различных по своим ФМХ кон- 
тактирующих пар как сталь-медь и сталь-полиэти- 
лентерефталат. В ячейках у контактной поверхности 
величины осевых и радиальных деформаций пары 
сталь-медь одного порядка. Для пары сталь-полиэ- 



тилентерефталат осевые и радиальные деформации 
трубы МУ становятся на порядок меньше, нежели в 
случае медного ВУ. Поэтому можно считать, что зна- 
чительного (по сравнению с размерами ячейки ко- 
нечно-разностной сетки) расхождения узлов двух тел 
в осевом направлении на поверхности контакта нет. 



4,3 4,4 4,5 4,6 4,7 4,8 4,9 5,0 

Рис. 6. Распределение реакции ведущего устройства по его 
ширине для: 1) медного ВУ при отсутствии снаряже- 
ния; 2) медного ВУ; 3) ВУ из полиэтилентерефталата; 
4) ВУ из ортотропного материала 

Деформации и напряжения в трубе МУ и МЭ 
для граничных условий (6) имеют несколько мень- 
шие значения. При таком способе закрепления 
торцов материал снаряжения начинает больше 
воспринимать нагрузку, несколько разгружая тем 
самым корпус МЭ под ВУ. В дополнение к зонам 
пластического деформирования, имеющим место 
для случая задания граничных условий (5), пласти- 
ческие деформации возникают также в снаряже- 
нии у оси вращения для і~Ь/2. Материал же трубы 
МУ для обоих случаев закрепления торцов не полу- 
чает пластических деформаций, что соответствует 
требованиям эксплуатации. Характер распределе- 
ния искомой реакции р(і) на контактной поверх- 
ности качественно остается таким же, как и при 
жесткой заделке торцов, но количественно ее вели- 
чина становится несколько меньше. Так, средняя 
величина реакции ВУ / р «^492 МПа. Таким обра- 
зом, использование граничных условий (5) и (6) 
позволяет получить своего рода крайние значения 
для оценки реальных параметров НДС. 

Расчет реакции р для этой модельной задачи по ме- 
тодике А.А. Ильюшина дает величину р т =-А 75,6 МПа, 
которая меньше аналогичных величин, полученных с 
помощью численной методики для условий закрепле- 
ния торцов (5) и (6) на 12,0 % и 3,5 % соответственно. 
Реакция р т вычисляется с использованием уравнений 
теории оболочек, гипотезы которой приводят к ужес- 
точению конструкции и, таким образом, к некоторому 
занижению величин определяемых параметров НДС. 
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Описанные численные результаты получены с 
погрешностями е=2,5 % выполнения равенства (2) 
и е„=5 % для вычисления перемещений в двух 
последовательных упругопластических задачах. 
Для получения решения задачи с заданной точ- 
ностью корректировка 21 радиальных перемеще- 
ний щ, щ на контактной поверхности "труба МУ - 
ВУ" по формулам (3) и (4) проводилась 14 раз. О 
степени выполнения статических граничных усло- 
вий можно судить на примере радиальных напря- 
жений ст г Следует отметить, что одним из досто- 
инств ВРМ решения задачи определения НДС яв- 
ляется автоматическое выполнение статических 
граничных условий (также как и выполнение урав- 
нений равновесия) при минимизации функциона- 
ла полной потенциальной энергии системы. 

4. Выводы 

Построенная достаточно строгая и полная в ма- 
тематическом смысле модель и созданная числен- 
ная методика позволяют определять реакцию веду- 
щего устройства на основе информации о физико- 
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механических характеристиках материалов и гео- 
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